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Progress in spintronics has been aided by char- 
acterization tools tailored to certain archetypical 
materials [1—6]. New device structures and mate- 
rials [7—12] will require characterization tools that 
are material independent, provide sufficient res- 
olution to image locally-varying spin properties 
and enable subsurface imaging. Here we report 
the demonstration of a novel spin-microscopy tool 
based on the variation of a global spin-precession 
signal in response to the localized magnetic field 
of a scanned probe. We map the local spin den- 
sity in optically pumped GaAs from this spatially- 
averaged signal with a resolution of 5.5 pm. This 
methodology is also applicable to other spin prop- 
erties and its resolution can be improved. It can 
extend spin microscopy to device structures not 
accessible by other techniques, such as buried in- 
terfaces and non-optically active materials, due 
to the universal nature of magnetic interactions 
between the spins and the probe. 

Precession in a magnetic field is a hallmark charac- 
teristic of a magnetic moment, which makes it a power- 
ful discriminant of spin-related phenomena in spintronic 
studies; this is exploited in Hanle effect measurements 
[13]. It can also result in loss of spin information in 
the presence of unwanted magnetic fields, especially spa- 
tially varying ones. However, in this report we show that 
the spin precession due to the spatially inhomogeneous 
magnetic field of a micromagnetic probe (pP) can be em- 
ployed to encode local spin information into a spatially 
integrated measurement. This information can then be 
decoded using standard deconvolution techniques. 

To demonstrate our technique we first measure a sig- 
nal proportional to the globally- averaged spin density, E, 
as a function of the /iP's position. From this measure- 
ment we can then obtain a quantity proportional to the 
unperturbed spin density (i.e. the steady state density in 
the absence of the probe), p, by deconvolving the signal 
from a theoretical or experimentally-determined Preces- 
sional Response Function (PRF). The PRF captures the 
response of the spins to the magnetic environment they 
experience. We repeat this process in the presence of 
an applied uniform transverse magnetic field to further 
demonstrate our understanding of this microscopy tool. 

We generate p = p(r s )z, in a GaAs membrane via 
optical pumping [13], where r s = (x s ,y s ,0) is the po- 
sition within the sample which we treat as two dimen- 



sional. Due to precession the steady state spin density 
in the presence of any magnetic field will be different 
from p, and we will denote it by S. We measure a glob- 
ally averaged spin-PhotoLuminescence (PL) signal [13], 
^ S z d 2 r s , where S z refers to the z component of 

S. 

This signal is measured while the pump laser is scanned 
relative to the \iP in the x and y directions. For our 
uniform sample, this is equivalent to scanning the pP 
relative to a fixed pump, and henceforth we will regard 
this to be the case. The field B p of the pP modifies the 
precession behavior of the spins (see Fig. 1). This field 
depends on /iP's position. A spatially uniform trans- 
verse magnetic field, Ptx, may also be applied to further 
tailor this precession. More details of the experimental 
set-up and measurement techniques are presented in the 
Supplementary Information (SI). 

Fig. 2b shows the measured PL intensity (which is pro- 
portional to p) for a particular p = p c , which will be used 
later for obtaining the PRF. Panels c and d show the cor- 
responding E c , for B t = 0T and 0.145 T respectively, as 
a function of the /iP's position, r p = (x p , y p: z p ). 

An expression for E in the limit of small diffu- 
sion (which our data shows is a reasonable approximation 
for this experiment; more details of the derivation in the 
SI) is given by: 

POO 

E(r p ,P t )oc / ^(r s ,r p )d 2 r s 

J oo 

rOO 

<x / H B (R,B t )p(r s )d 2 r s 

J OO 

= H B (R,B t )*p(r s ) (1) 
where, * represents a convolution, R = r„ — r s and 
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Hb is the PRF and is the signal which would be col- 
lected for a spatial delta function injection. B± and Bn 
are respectively the magnitudes of the perpendicular and 
parallel components of the total field, B = B p (R) + P^x, 
experienced by the spins. The components are defined 
with respect to z, the orientation of injected spins. The 
gyromagnetic ratio is denoted by 7 and r s is the spin 
relaxation time. 

As seen from the previous equations, S z is decreased 
by B± because it causes the spins to precess away from 
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FIG. 1. The simulations shown illustrate the key physics underlying Spin Precession Imaging; the inhomogeneous field of 
a micromagnetic probe pP (sphere with black arrow, indicating the magnetization direction) generates a well-understood 
spatially- varying dephasing of spins. Spins that are injected into a semiconducting sample reach a steady-state density resulting 
from a combination of the local spin properties, the injected density, and the magnetic field due to the pP. The arrows in all 
figures represent the steady-state spin density vectors S. The color scale represents the z component (parallel to the orientation 
of the injected spin) of the spins, a) Demonstrates the case of a spatially-uniform injection density and highlights the full 
spatial-extent of the /xP's influence. Strong transverse fields from the probe dephase spins in an annular disk centered beneath 
the probe. Parallel fields immediately below the probe protect spins from dephasing. b) and c) show simulations identical to 
panel a) except for the case of a spatially confined injection profile, corresponding to two different pP positions. Panels b and 
c correspond to the blue and red dots in Fig. 2 c respectively. 



the injected direction, resulting in a dephasing of the 
ensemble. On the other hand Bu keeps the spins from 
tipping away from the injection direction, resulting in a 
small Ob and large S z . We can view as an effective 
dephasing factor [14, 15]. 

B\\ and B± have distinct spatial variation in our exper- 
iments, and the consequences of their competing effects 
are evident in Fig. 3, where we show line scans (along 
x and y) for several values of B t . Also shown are fits 
obtained from Eqn. 1 in which the probe is modeled as 
a point dipole with a moment m = m p z located a height 
z p above the sample (see Fig. 3 caption and the SI for 
more details). 

The peak, marked by vertical green dashed lines at x or 
y = 0, occurs when the pP is located directly above the 
point of maximum injected spin density. At this point 
there is a maximum in Bn (~0.8 T) from the pP that 
preserves S z . When the pump is far from the /iP, the 
signal decreases with increasing B t with a Lorentzian line 
shape (as expected in a Hanle measurement) whose half- 
width, B1/2 = (t^s) -1 , is 0.0111 T for our experiments. 

A second peak (blue dashed line) seen in the line scans 
along x occurs where the x component of the field from 
the /iP cancels B t . As B t is increased, this point occurs 
closer to the \iP where its field is stronger. 

The fits indicate the effectiveness of Eqn. 1 in describ- 
ing our data. The sensitivity of the global signal to spins 
at different locations relative to the /iP, described by the 
convolution, forms the basis for imaging. 

To obtain an unknown spin density, p Ul from our mea- 
sured signal, the PRF needs to be known. The PRF can 
be obtained theoretically or from experimental data, if 
we have a known p. For the latter, we use the camera 



data shown in Fig 2d as being proportional to p c . Then 
ife(R, Bt) — S c (r p , B t ) ® p c (r s ), where we use ® to in- 
dicate a deconvolution process. We use the Wiener algo- 
rithm [16] to implement the deconvolution. The resulting 
experimental PRFs for both low and high B t are shown 
in Fig 4. For comparison the theoretical PRFs are also 
shown; more details for obtaining them are presented in 
the SI. 

To test the fidelity of our imaging process, we now 
use the experimental PRF for an unknown p u . The 
measured S n , at low and high fields, is presented in 
panels a and c of Fig. 5. We then extract p u (r s ) — 
£ n (r p ,P t ) ® iTe(R, Bt). The extracted spin densities 
are shown in panels b and d of the same figure. Also 
shown (panel e), for independent verification of our imag- 
ing technique, is a camera image for the PL (oc p u ). The 
linecuts present a more quantitative comparison of the 
extracted and measured data. 

The ability to extract the spin density with both high 
and low B t shows the exclusion of spurious effects, such 
as reflectivity changes as a function of the /iP's position, 
in our data. Also, B t provides a knob to optimize the 
PRF to suit particular imaging needs. High field imag- 
ing might provide a more intuitive PRF for the case of 
global detection. Low-field imaging maybe more useful 
for non-local electrical devices, where a large transverse 
field would dephase spins before they reach the detector. 

The line scans shown in Fig. 3 provide us a measure of 
the spatial resolution, £, in our experiment. A Gaussian 
fitting of the narrowest lobe gives us £ < 5.5 pm. It 
is an upper bound that is being set by the feature size 
we are imaging. As in magnetic resonance imaging, the 
magnetic field gradient, ft, sets the ultimate resolution in 
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FIG. 2. Spatial maps of the injection beam profile from the optical pump, the injected spin profile and the response of the 
global spin signal to probe position relative to the pump, a) Spin-insensitive camera image showing the PL spot created by the 
pump beam in a 1 /im thick GaAs membrane. Also visible is the the NdFeB micromagnetic probe fiP, magnetized along the 
sample normal and glued to the membrane on the side opposite to that of the pump beam. An additional uniform magnetic 
field Bt x may be applied. In the experiment the pump beam is scanned relative to the jiP using an objective mounted on 
translation stages, b) Camera image of PL that corresponds to (and is proportional to) injection profile p c . c) The measured 
spin signal S c (colorbar) corresponding to the injection profile from panel b, at B t — 0. The location of each pixel corresponds 
to the relative position between the jiP and the pump beam; a blue and red dot are included to provide two example cases. 
The location of the blue dot corresponds to a pP position directly above the injection beam; this configuration results in a spin 
density mostly pointing along the z direction (Fig. lb) and gives a large signal. The red dot corresponds to a fiP position 8 
/xm away from the center of the beam; spins precess away from the z direction due to large perpendicular fields (Fig. lc) and 
thereby reduce the signal, d) E c for a large B t = 0.145 T. 



the absence of diffusion, ( = B 1 / 2 /k [14]. Gradients of 
up to ~ 4 x 10 6 T/m have been reported recently [17], 
which are at least an order of magnitude larger than in 
this experiment, and should enable much finer resolution. 

Resolution will be limited by diffusion and the un- 
avoidable reduction of the spin signal as the detected vol- 
ume shrinks. While diffusion can degrade the resolution, 
the spatial precessional response can be numerically an- 
alyzed to obtain valid and useful spatially resolved data. 
For large enough gradients, sub-diffusion length and sub- 
diffraction limit resolution should be achievable. Images 
are obtained in the presence of spin diffusion by MRI [18- 
20] ; this should be feasible for spin precession imaging as 
well. 

In summary, we have demonstrated a new technique for 
imaging spin properties using the precessional response 
of spins to a micromagnetic probe's field. While we have 
imaged the variations in the spin density, the technique 



is more general since the response of the spins is sensitive 
to a variety of spin characteristics including spin lifetime 
and gyromagnetic ratio. Work is underway to general- 
ize the technique presented here using scannable probes 
mounted on cantilevers. Due to the magnetic nature of 
interaction between the probe and the spins, which can 
extend through layers of a heterostructure and a few mi- 
crons deep, this tool should enable subsurface imaging. 
This technique should be applicable to a wide variety of 
materials because it relies on proven spin polarization 
detection techniques. With optical detection it can en- 
hance imaging resolution, and with electrical detection it 
can enable imaging where none exists at present. 

Funding for this research was provided by the Center 
for Emergent Materials at the Ohio State University, an 
NSF MRSEC (Award Number DMR-0820414). We wish 
to thank Cristian Cernov for creating the rendered im- 
ages presented in this article and the SI. 
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FIG. 3. Spin signal (E c , corresponding to p c ) for various B t plotted as a function of the relative position between the pump and 
the pP along the y (left) and the x (right) directions. The data at various Bt are offset vertically for clarity of presentation. 
The open circles are experimental data while the solid lines represent fits obtained from Eqn. 1 using m z — 2 x 10 -9 J/T, 
z v — 8 /im, and r s = 2.33 ns. The spin lifetime r s is obtained from Hanle measurements (see SI) and the pP radius (z p ) is 
taken to be 9 /xm. The peaks marked by the green dashed lines result when the peak in p c lies directly under the pP. The net 
parallel field directly under the fiP is mostly due to B pz (the z component of the /xP's field) and exceeds the net transverse 
field (mostly given by Bt); this results in a large E (Eqn. 1). The second set of peaks seen in scans along the x direction (blue 



dashed line) occur when B p 



-B t (B px is the x component of /iP's field). 
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FIG. 4. a) Precessional Response Function obtained through Wiener deconvolution of p c (Fig. 2b) from E c (i^ = 0) (Fig. 2c). 
b) Theoretically derived PRF obtained from Eqn. 1, using the dipole moment from the fits in Fig. 3. c) and d) similar 
experimental and theoretical PRFs for B t — 0.145T. The experimental PRF data has been normalized and offset to highlight 
the match between theory and experiment. 
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FIG. 5. Obtaining the spatial variation of an unknown spin density from the measured signal E: a) Spin signal S w measured 
for a spin density profile, p u (more information on how it was produced given in the SI), with Bt — 0T. b) An image of p u 
extracted from the Wiener deconvolution of Y, u (panel a) using the experimental HsiBt = 0) (Fig. 4a) as the deconvolution 
kernel, c) and d) Similarly measured H u and extracted p u respectively for B t — 0.145 T. e) Line cuts of p u taken along the 
dashed lines in panels b (blue) and d (green) . Also shown in red is a line cut from an independent camera image of p u shown 
in panel f. 
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I. SAMPLE 

The sample is a 1 jam thick (and ~ 3 mmx ~ 3 mm) n-GaAs (001) membrane prepared by etching away the substrate 
underneath an MOCVD grown epitaxial layer. Fig. SI shows the growth structure, where the top layer forms the 
membrane. A nominally lattice-matched InGaP layer is grown under this device layer to act as a stop-etch layer 
during the fabrication process. SIMS analysis (EAG labs) shows that the doping level of Si in the device layer is 
~ 1.4 x 10 22 m -3 . This is below the metal-insulator transition level for GaAs. Thus the samples are insulating during 
our low temperature (< 20 K) measurements. This supports the lack of significant diffusion in our measurements. 

To prepare for the etching, the device layer side was glued to a 0.5 mm thick single crystal (OOOl)-Sapphire sub- 
strate (5 mm x 5 mm), whose c-axis orientation was chosen to reduce birefringence during the polarization dependent 
measurements. The sample was glued using Epotek 301-2 optical epoxy to reduce PL background. 

The substrate was etched away following a recipe given elsewhere 1 . The etch stops at the InGaP layer because of the 
high etch selectivity. The InGaP is then etched away using HC1. After the etching process, a NdFeB micromagnetic 
particle was glued to the membrane (on the side opposite to the sapphire and where the substrate used to be) using 
micro-manipulators, under an optical microscope. 




GaAs semi-insulating substrate 



I I 

FIG. SI. The MOCVD growth structure for the GaAs sample. 



II. THE OPTICAL SET-UP 

Spins are injected into the GaAs membrane using optical spin injection, which is a well established technique. 
The detection is based on measuring the circular polarization of the PL, which in turn is proportional to the spin 
polarization within the sample. The set-up used for optical pumping and the spin-PL measurements is shown in 
Fig. S2. Also, provided in Table SI is a detailed list of the various optical components used and their purpose. 

The Electro-Optic Modulator (EOM) allows us to modulate the polarization of light, and thus modulate the spin 
polarization within the sample, between +z and — z, and enable lock-in measurements. The EOM modulates the 
780 nm pump between two linear polarization states. These are converted into the two circular states by the the 
Quarter Wave Plate (QWP). The circularly polarized laser light is focused on the sample (which sits inside an optical 
cryostat at a temperature of 17 K) using an objective. A wire mesh is placed in the path of the pump beam to 
produce the spin density profile p u (see main text). When producing different spin density profiles, the pump power 
was adjusted to keep the Hanle half- width far away from the jiP equal. 

The circularly polarized light injects spins due to a combination of spin-orbit coupling and optical selection rules. 
The resulting PL due to the recombination of spin-polarized carriers is circularly polarzied. This PL (at the band 
edge of GaAs, 819 nm) is then collimated by the objective and converted back into the linear states by the QWP. 
The Wollaston prism splits the two orthogonal linear states; each of which is then collected by a separate photodiode 
that are part of a diode-bridge circuit. More details of how the photodiode signal is processed is presented in the next 
section. 

The end of a multi-mode fiber is placed ~1 mm away from the sample (on the side opposite to laser injection, i.e. 
on the \iP side). The fiber is used for illuminating the sample with unpolarized broadband light (from a halogen 
lamp) to assist with camera imaging and tracking of the fiP during the measurements. The back illumination gives 
negligible spin signal for our lock-in measurements. 
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S.no 


Description 


Comments 


A 


GaAs sample 


see section I for more details 


B 


Electromagnet 


Power supply controlled via DAQ card 


C 


Optical crystat 


Operated at 17K 


D 


Objective lens 


plan fluor 10X 


E 


Quarter-wave plate 




F 


Dichroic mirror 


Reflects pump towards the sample, but attentuates it during transmission 


G 


Notch Filter 


For 780 nm 


H 


Cube beam splitter 


Non-polarizing 


I 


Focusing lens 


f 100 mm 


J 


Wollaston prism 


Separates the two linear polarization components 


K 


Bandpass filter 


Centered at 820 nm, 10 nm width 


L 


Photodiode bridge 


The difference channel is amplified by a factor of 100000 


M 


focusing lens 


f 200 mm 


N 


Cube beam splitter 







Long pass filter 


cut off at 850nm 


P 


CMOS camera 




Q 


Photodiode 


To measure spin-insensitive DC level of the PL 


R 


Wire mesh 




S 


Electro-optic modulator 


Modulated polarization of light 


T 


Polarizer 


transmits vertical polarization, as needed for the EOM operation 


U 


Polarizer 


Allows for fine control of pump power 


V 


ND filters 




W 


Pump laser 


780 nm, stable laser source, intensity modulated 


X 


Function generator 


used to provide modulation signal for EOM, square, fmod ~ 73 Hz 


Y 


Voltage preamplifier 


Amplifies the A-B channel of the dide bridge, Gain of 200 


Z 


Voltage preamplifier 


Amplifies the B channel of the dide bridge, Gain of -1000 


AA 


Voltage preamplifier 


Amplifies the A channel of the diode bridge, Gain of 1000 


AB 


Lock-in amplifier 


Demodulates the pump power modulation, typical time const. 200-500ms 


AC 


Lock-in amplifier 


Demodulates the EOM modulation, typical time const. 200-500ms 


AD 


Function generator 


used to provide modulation signal for pump intensity, square, fmod ~ 1100 Hz 


AE 


Computer 


Controls the instruments and performs data acq., thru GPIB and DAQ card 



TABLE SI. A list of the various optical components and electronic instruments used in the experiment. The s.no correspond 
to the letters in Fig. S2 and S3. 

Also, a variety of optical filters are used in various parts of the optical set-up to remove unwanted light from reaching 
the photodiodes or the camera. Please see Table SI for more details. 



III. SCANNING AND DATA COLLECTION 

The data shown in Fig. 2b, c and Fig. 4a, c of the main text are obtained by scanning the objective (which in turn 
scans the pump beam relative to the sample) using motorized translation stages. The position of the jaP relative to 
the beam was tracked using a home built software solution. The Labview-based software uses the camera image and 
pattern recognition algorithms from National Instruments. This measured position in pixels of the camera image was 
converted into microns using a SEM image of the \iP as a calibration. 

Fig. S3 shows the instrumentation and hardware signal processing used for obtaining the spin signal from the 
photodiode voltages. The signals from the photo-diode bride allow us to compute the average steady state spin 
polarization within the sample, £ ex y^qr^S where the subscripts refer to right or left circular polarization of the PL. 
The difference signal, Vr — Vl, is the lockin signal due to the modulation of the polarization using the EOM. Vr + Vl 
is measured as a different lock-in signal by modulating the power of the pump laser. The normalization by the sum 
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FIG. S2. The optical set-up for the optical pumping and spin PL detection. Please see Table SI and sec. II for more details of 
the various components. 



of the two circular components removes spurious reflectivity changes in the signal, as seen in Fig. S4. 

The step size of the translation stages is somewhat variable. To allow for further data analysis, we have interpolated 
the 2D scan data (but not the line scan data) using Igor Pro software to have the signal on a uniform spatial grid. 
The data shown in the main text is the interpolated data. Further, all the data has been normalized such that E 
far away from fiP at zero field and maximum field are taken to be 1 and respectively. Also, for the line scans the 
position values were offset to bring the central peak to xoiy = 0. 



IV. SPIN DYNAMICS IN A NON-MAGNETIC MATERIAL 

In a two-dimensional non-magnetic semiconductor, the spin density may be governed by the following equation of 
motion 2 

<9S S 

— = G + DV 2 S + C(E • V)S + 7 B x S , (SI) 

Ot T s 

where all vector fields (S,G,E,B) are, in general, functions of the spatial coordinate r s = (x s ,y 3 ,Q) within the 
sample and time t. The spin density is given by S; each component of this vector field gives the difference between 
spin "up" and "down" particle densities in that direction. The first term on the R.H.S. of Eqn. SI, G, represents 
the rate at which spin density is externally injected and related to the intensity of the pump light. The second and 
third terms dictate the diffusion and drift of the spins respectively, where D is the diffusion constant, E is the electric 
field, and f the electron (or hole) mobility. We assume no difference in diffusivity between opposite-spin carriers in 
any direction. The fourth term gives the precession of spin density around a net magnetic field B, with 7 = g/ie/Zz 
being the gyromagnetic ratio; where g is the g-factor, \ib is the Bohr magneton, and h the reduced Planck constant. 
The final term represents spin relaxation with a spin lifetime of r 3 . 

We study the case in which E = 0, r s is spatially uniform, and B and G vary spatially but not in time. We also 
consider an injection rate constrained to the z-axis, i.e. G = G z z. From a device point-of-view, one usually measures 
the steady-state solution of spin density S(x,y,t — >• 00); in what follows we will only consider this solution. 

In the limit where diffusion is negligible, Eqn. SI possesses a purely algebraic solution, which can be written more 
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FIG. S3. Schematic of the signal processing done to obtain the measured signal. Please see Table SI and sec. Ill for more 
details of the various components. 



conveniently by introducing the following notation 



/ B* -B* \ 



B 



-B* z 



B* 
-B* / 



where we have applied a scaling to B (and its vector components B x ,B y ,B z \ B* = 7T S B . 
The steady state solution, dS/dt = 0, for the case of D = is then given by 

S(r) = [/-B(r,)]-V(r s ) 



(S2) 



where, p(r s ) = G z (r s )r s . It should be noted that p is the steady state spin density in the absence of any magnetic 
field and is the quantity that we are interested in measuring in this experiment from a signal that is averaged in r. 

Seen in the above equation is that fact that the absence of both drift and diffusion eliminates all coupling between 
neighboring positions (which in general are coupled due to derivative operators). Thus, in this particular limit, the 
steady-state solution is determined algebraically by the magnetic field values at a given position, and does not need 
to be solved using the full differential equation. The solutions, based on Eqn. S2, for the various vector components 
of S(r s ) are determined by B and p only at that r s ; they are given by, 



Sy — 



BIBt - B* 



1 + Bf + Bf + 




y z + B x 




1 + B* 2 + Bf + 


Bf" 


1 + Bf 

"1 i 7~> 2 i 7~~> -J-- 2 i 





(S3a) 
(S3b) 
(S3c) 
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V. THE MEASURED SIGNAL 



In our experiment the measured signal may be described by 

Sex / S z {v)dx s dy s , (S4) 

J A 

where A represents the area of detection and can be assumed to extend to infinity. The magnetic field, B, in our 
experiment is provided by the field from a micromagnetic probe (/iP), B p , and a spatially uniform field PtX. Assuming 
the \iP to be a point dipole, the field experienced by the spins, at position r s within the sample, is given by 



B(R) = B P (R) + B t x 

/i 3R(m • R) - mi? 2 

4^ W> 



• B t x. (S5) 



Here, R = r s — r p is the relative spatial vector between the position of the /iP, r p , and the position of the spins within 
the sample; m is the moment of the jiP and /xq is the permeability of free space. 
Substituting Eqns. S3 and S5 into Eqn. S4, we can write the measured signal as 

POO POO 

£(r p ) oc / / S z (Y Sl Y p )dx s dy s 

J oo J oo 

X L L {l + B* A Kf + B;(Kf + B*(Iif) p{rs)dXsdys 

<xH B (R)*p(r s ). (S6) 

The ★ operator denotes a two-dimensional convolution integral and H B is the Precessional Response Function (PRF), 
which may be written in the following manner to emphasize the physics at play 

H (R) = ± + B?(R) 

= 1 + 0|(R)' (ST) 
where the effective dephasing factor 9b is given by 

n2 (n) _ l 2 B±m 2 r! 

^ (R) - 1 + 7 2B || (R)2 T |' ( S8 ) 

where B\\ = yjB% + P^ and B± = B z refer to the parallel and perpendicular components of the total field B, with 
respect to the injected spin direction z. 



VI. FITTING THE LINE SCANS 



The fits to the line scans shown in Fig. 3 of the main text were done using the convolution equation shown in 
Eqn. 1 of the main text (same as Eqn. S6). The various parameters for the fit were either measured independently or 
constrained by measured quantities. The P1/2 = l/7 r s was obtained using Hanle measurement done far away from 
the jiP. This Hanle data is shown in Fig. S5, along with the Lorenztian fitting done to obtain the halfwidth. The 
fitted P1/2 = 111G corresponds to a lifetime of 2.33 ns assuming a g-factor of -0.44 for GaAs. No diffusion (D = 0) 
was assumed. A moment of 2 x 10 -9 J/T and a moment height of 8 fim were used. These parameters provide a better 
fit (in terms of mean squared error) than 3 x 10 -9 J/T and 9 /im, which would have been nominally expected values 
given the particle radius (9 /im) as determined from SEM images. This may be due to the limitations of the simple 
single dipole model that we use for fiP. The injection spot is assumed to be a Gaussian of 5 fim half- width, which is 
close to the the values obtained from fitting the camera image data to a 2D Gaussian. 
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FIG. S4. Background removal and normalization of spin signal E. All the data shown in the figure correspond to the line 
scan data shown in Fig. 3 of the main text (Bt — 0, scan along x). Two lock-in signals A — B (blue dashed line, right axis) 
and A+B (green dashed line, right axis) are measured, where A and B correspond to the two photodiodes in the diode-bridge 
circuit. These correspond to the difference and sum of the intensities of the two circularly polarized components of the PL. The 
spin density is then proportional to A-B/A+B (red solid line, left axis). The A-B and A+B channels can exhibit reflectivity 
changes (for e.g., the dip at yp 20 /mi) as a function of the probe's position or due to the surface conditions of the sample. 
However, these are reasonably removed during the division process. Further evidence for the spurious reflectivity changes is 
seen in the spin-insensitive PL signal (black dashed line, right axis) collected by a separate photodiode. 
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FIG. S5. Hanle curves measured far away from the pP for the two spin density profiles p c (red circles) and p u (blue pluses) 



VII. EXTRACTING THE SPIN DENSITY 



The PRFs presented in Fig. 4a and b were obtained by a Wiener deconvolution process, 

H B (R,B t ) = X c (r p ,B t )®p(r c ). 



(S9) 



The deconvolution was implemented using Mathematica with a regularization parameter of 20. 

The PRFs presented in Fig. 4c and d were calculated using Eqns. S7, S8 and S5, and the point dipole parameters 
used in sec. VI. Fig. S6 presents the p u obtained by deconvolving using the theoretical PRFs in the main text from 
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FIG. S6. spin density p u obtained from deconvolution of T, u and the theoretical PRFs shown in Fig. 3a (left hand panel) and 
Fig. 3b (right hand panel) of the main text. 



p(r c ) = E c (r p ,B t )®H B (R,B t ). 
These deconvolutions used a regularization parameter of 10. 



(S10) 
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